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Abstract 

We prove that the number of 3-dimensional simphcial complexes having the 
spherical topology grows exponentially as a function of a volume. It is sug- 
gested that the 3d simplicial quantum gravity has qualitatively the same 
phase structure as the 0(n) matrix model in the dense phase. 



Present address. 



1 Introduction 

The success of the matrix models as a theory of 2D quantum gravity and non- 
critical strings has brought about a hope that analogous discrete approaches 
might be instructive in higher dimensions as well. The most natural way 
to introduce discretized quantum gravity is to consider simplicial complexes 
instead of continuous manifolds. Then, a path integral over metrics (which 
is the crux of any approach to quantum gravity) can be simply defined as 
a sum over all complexes having some fixed properties. For example, it is 
natural to restrict their topology. In the present paper, we consider only the 
3-dimensional case, where, by definition, a simplicial complex is a collection 
of tetrahedra glued by their faces in such a way that any two of them can 
have at most one triangle in common. 

To introduce metric properties, one can assume that complexes represent 
piece- wise linear manifolds glued of equilateral tetrahedra . The volume is 
proportional to the number of them. In a continuum limit, as this number 
grows, the length of links, a, simultaneously tends to zero: a — 0. In 
three dimensions, there is the one-to-one correspondence between topological 
classes of piece- wise linear and smooth manifolds, therefore the model is self- 
consistent. 

Within a piece-wise linear approximation, the curvature is singular: the 
space is fiat everywhere off links. Therefore, one should consider integrated 
quantities, the main of which is the mean scalar curvature 



a is the lattice spacing; N^, Ni, N2, are the numbers of vertices, links, 
triangles and tetrahedra, respectively. 

For manifolds in 3 dimensions, the Euler character vanishes x = ^0 ~ 
Ni + N2 — = 0. Together with the other constraint N2 = it means 
that only 2 of A^j's are independent, A'^i = Nq + N^, and a natural lattice 
action depends on 2 dimensionless parameters. Now, we are in a position to 
define the partition function for simplicial gravity p|, ^, ^ Q: 
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where, for concreteness, we have restricted the topology of complexes: Z]{S3} 
denotes the sum over all simplicial 3-spheres, S^. 

It can be shown that, for fixed, Nq is restricted from above as Nq < 
^N3 + const. Therefore, it is natural, keeping a fixed, to tend /i to its critical 
value, fic, when the sum over in Eq. (H) becomes divergent and the mean 
volume, {N^), tends to the infinity. In the vicinity of one could expect to 
find critical behavior corresponding to a continuum limit of the model. 

For this scenario, it is crucial that /ic is finite: < /ic < +00, i.e., the 
number of complexes as a function of the volume, N3, should grow at most 
exponentially The purpose of the present paper is to show that it is 
indeed the case. More precisely, we prove 

THEOREM: There exists a finite constant < /i* < +00 such that the number 
of 3- dimensional simply- connected simplicial complexes, AflN^), constructed 
of the given number of 3-simplexes, N^, obeys the restriction 

AfiNs) < e'^*^^ (3) 

for big enough values of N3. 

"Simply-connected" means that all loops in a complex C are contractible, 
i.e., Tti{C) = 0. If the Poincare hypothesis is true, this class of complexes 
coincides with homeomorphic spheres. If not, it is somehow wider. However, 
the statement of the theorem will obviously hold for the spheres as well. 

2 Proof of the theorem 

We want to estimate //{N^) from above being unable to calculate it precisely. 
The general principle can be formulated as follows. Let S{a) be a finite set 
of objects defined by a set of descriptions a. If a C (3, then S{[3) C S{a) 
and \S{j3)\ < |>S'(a)|, where \S\ denotes the number of objects in a set S. In 
other words, being less specific, one cannot decrease the number of objects. 

Throughout the paper we use the following trick: if S = U{a} ^a, then 
\S\ < |A|max{^}|i?^| and, if all Ba 7^ 0, then \A\ < \S\. Less formally, let us 
notice that the existence of the exponential bound (|) is not spoiled, if one 
weights every complex with a non-zero weight which itself depends on 
and grows at most exponentially. Also, if one manages to represent a sum 
over complexes as a weighted sum over another class of objects and to prove 
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that the weights are exponentially bounded, then one can put them all equal 
to 1 and simply estimate the number of objects in the new class. 

The simplest example of a weighted sum over complexes is given by the 
micro-canonical partition function Z^^la) in Eq. (^. Another example of 
physical interest is the partition function in the presence of free matter fields. 
To introduce it let us consider the set of cell complexes dual to simplicial ones. 
Their 1-skeletons are some <p^ Feynman graphs whose vertices correspond to 
tetrahedra and propagators are dual to triangles. Such a graph can be defined 
by the adjacency matrix 

^ _ f 1 if vertices i and j are connected by a link 
^■^ 1 otherwise 

Of course, the set of 0^ graphs is not identical to the one of simplicial 
complexes, similarly as, in two dimensions, the set of ordinary (p^ diagrams 
is different from the one of triangulations (one has to introduce the notion 
of fat graphs to establish the equivalence). However, we shall loosely use the 
term "0^-diagram" because of its commonness. 

If one attaches the n-component free matter field (/x = 1, . . . , n) to 
the i-th vertex, the resulting gaussian integral can be performed explicitly 
for every given graph G: 



/N3 — I n 1 iV3 — i 

n n dxf exp [ - - ^ G., (xf - x^y] = det "^L (5) 
i=i ^=1 ^ ij=i 

where L is the discrete Laplacian 

Lij = ASij - Gij (6) 

for the determinant of which there is the nice combinatorial representation 
given by the Kirchhoff theorem: 

detL=|r(G')| (7) 

where T{G) is the set of all maximal connected trees embedded into the (p^ 
graph G or, equivalently, the number of ways to cut links of the graph such 
that it becomes a connected tree. To kill the zero mode in Eq. (^, we fixed 
the field at the iVs-th vertex, therefore we consider rooted trees with the root 
attached to this vertex. 
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The case we are interested in is the 2-component Grassmann field, where 
we obtain the partition function^ 

{a, //) = E det L e"^^-'^^^ = J2 (a)e-^^^ (8) 

{S3} N3 

The weights, det L, are non-zero integers, hence, 

ZM,{a) < Z^^lia) (9) 

In two dimensions, this type of matter has the central charge c = —2 
and the corresponding matrix model has been solved explicitly We can 
repeat the same trick in 3 dimensions: 

4i(«) = EdetLe"^«=E e"^« = J] Yl (W) 

{G} {G} {T(G)} m {G(T)} 

Namely, we have the sum over complexes (0^ graphs), J2{G}y weighted with 
the number of maximal trees in each J2{t{g)}- We can obtain exactly the 
same configurations taking the sum over all possible trees, J2{T}y weighted 
with the number of complexes, J2{g{t)}j which can be recovered from a given 
tree, T, by restoring cut links. 

In 3 dimensions, the appearing trees are dual to spherical simplicial balls 
obtained from a single tetrahedron by subsequently gluing other tetrahedra 
to faces of the boundary. If we denote no, ni, n2 the numbers of vertices, 
links and triangles on the boundary of a ball, we find that 

no = iV3 + 3 n2 = 2(iV3 + l) (11) 

The set of the boundaries are a subset of all planar spherical triangula- 
tions, and their number is exponentially bounded as a function of A'3. There- 
fore, the remaining problem is to estimate the maximal number of closed 
3-dimensional simplicial spheres which can be obtained from an arbitrary 
spherical ball by pair-wise identifying triangles on its boundary. 

It seems to be rather difficult to find a criterion which would single out 
only homeomorphic spheres from all possible configurations. Therefore, we 

This model was considered by V.A.Kazakov several years ago as the simplest gener- 
alization of the corresponding 2-dimensional matrix model ||^ . 
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demand instead that all final complexes are to be simply-connected, i.e., 
having the trivial fundamental group. 

Let us consider loops embedded into a 0^ diagram dual to a simplicial 
complex C (which is the 1-skeleton of a dual cell complex, C). As any 
closed loop is contractible, it can be represented as the boundary of some 
2-dimensional disc embedded into the complex C. It means simply that, to 
shrink a loop, one has to pull it subsequently through 2-cells of C (they are 
in correspondence with links in C). Obviously, we have to allow for self- 
intersections and multiple coverings of the same face. However, any A;-cell 
of the disc has to cover one and only one fc-cell of C. In other words, we 
consider only regular {i.e., locally unambiguous) embeddings. 

If two triangles on the surface of a ball are identified, there is always a 
closed in C loop represented by a segment lying completely inside the ball 
with the ends attached to the identified triangles. Having cut C, one cuts 
any disc spanned by the loop. The discs are represented by planar graphs 
and the cut links form what is called "rainbow diagrams", i.e., collections 
of self-avoiding segments. If a cut disc remains connected, these segments 
form rainbows consisting of embedded arches. The deepest segment con- 
nects a pair of triangles having a common side on the surface of the ball, 
thus recovering a 2-cell of C. If the disc is disconnected by cutting C, such 
configurations may be absent. 

Proposition 1 Let C\, he a cell complex dual to a simplicial hall (i.e., oh- 
tained from some closed complex C by cutting its links). If, hy cutting a link 
in Cb, one disconnects all discs embedded into Cb and spanned by loops going 
through the link, while the complex remains connected, then C^, is not simply 
connected. 

Proof: If a cut disconnects all the discs, it means that the link enters at 
least twice in any closed contractible loop going through it. As the complex 
remains connected, there must be another path connecting the ends of the 
link. The union of this path and the link represents a closed loop which 
cannot be a boundary of any disc, hence, it is not contractible. 

Now, let us notice that, if one starts with an arbitrary closed simply- 
connected complex, a final ball does not depend on a sequence in which 
cuts are performed. One can choose a particular sequence never creating 
non-contractible loops, owing to the following fact. 
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Proposition 2 Let Cb be a cell complex defined above. If, by cutting any link 
in Cb, one creates a non-contractible loop, then Cb is not simply- connected. 
Proof: Let a cut create a new non-contractible loop. It is a collection of 
links of Cb- Before the cut, it was contractible, therefore there was a disc 
spanned by it. The cutting made it a ring, hence, the cut link belongs to 
the disc. Moreover, it belongs to any disc whose boundary is homotopic 
to the new loop. As it holds for any link, we are considering an effectively 
2-dimensional object, i.e., a collection of surfaces with holes having some 
identified vertices, links and/or faces. If, by deleting a link, we create a non- 
contractible loop in a simply-connected 2-dimensional space, it means that 
it has the homotopy type of a disc. A disc must have a boundary. However, 
we cannot create a non-contractible loop by cutting a link belonging to the 
boundary. Hence, we find a contradiction: one cannot fulfill all the conditions 
simultaneously. 

Thus, any simply-connected simplicial complex can be obtained from a 
ball by subsequently gluing up triangles having a common side on the sur- 
face and such a pair always exists. Obviously, the number of inequivalent 
sequences of such identifications is bigger than the actual number of com- 
plexes: two triangles may have two common sides and it may happen at any 
step. It means an overcounting. However, the gain is that we obtain a purely 
2-dimensional problem now. 

The boundary of the ball can become a collection of disconnected 2- 
dimensional spheres but any identification of triangles belonging to different 
components produces a non-contractible loop. 

So, one starts with an arbitrary 2-dimensional spherical (p^ diagram and 
applies the operation which can be represented as the flip of a link with the 
subsequent ehmination of it: 




Both of these moves are well known and have been used in Monte-Carlo sim- 
ulations of triangulated surfaces 0. Having glued all triangles, one flnishes 
with a collection of self-avoiding closed loops, the number of which equals 
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Let us notice that, instead of removing links, we can mark them with 
dashed propagators. For it we need the infinite set of vertices 



and so on 



(13) 



which are generated by flips. For example, 




and so on 



and so on. 

In the end, we obtain planar diagrams with two types of propagators: 
solid ones produce closed loops while dashed form arbitrary clusters. To 
recover initial diagrams one has to flip dashed links in all possible ways. For 
any given cluster, the number of flnal conflgurations is a subset of all planar 
(f)^ diagrams, hence, is exponentially bounded as a function of the number of 
links. Let us show it. At first, fiips have to be performed in linear clusters 
which gives all possible planar cjy^ trees doubled in a "mirror" : 
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and so on. The number of configurations with n dashed links, C„, is given 
by Catalan's numbers generated by the function 



n+l 



C{x) = Y.C^x^ ( T 



n=0 



n=0 




2x 



(14) 



After that dashed propagators attached to the ends of the hnear cluster 
have to be reattached in all possible ways to links of these "doubled trees" . 
The number of configurations can be estimated as follows. 

First, we have to calculate the generating function for the number of 
possible attachments of dashed propagators to a solid line: 




(15) 



It can be obtained from the equation 

q{iy,IJ.) 



+ zr^ 



1 — // 1 — u 



q{u, /x) + 1 



(16) 
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which gives 



(17) 



[1 — u — — (z/ + yu)z//i + v'^ij? 
Second, we have to know the generating function for the number of all 
possible attachments of dashed propagators to a /c-legged tree on the left (or 
on the right but not on both sides): 



fc+i 



t{x, v) n{ :;ii(T) ^(^' 0) = C{x) 



(18) 



The equation for it is simply 



which gives 



xC{x) 



Then, the quantity we need to calculate is 




(19) 



(20) 




n,m,k 



n,m,k 



k+1 




x{l-vY{l-iif 



((1 - v f - xC{x)) ((1 - - xC{ 



X. 



(21) 



Newly appearing clusters have to be expanded in all possible ways ac- 
cording to the same procedure. What is important is that their number is 
restricted by the coefficients of //), which is a factor in (PT]), and does 
not depend on a volume of the trees. Therefore, one simply iterates till all 
dashed links disappear. The combinatorial reason for the existence of an 
exponential bound for the number of obtained configurations is clear: they 
are all planar! 

To calculate the number of "unffipped" configurations, we can use the 
following matrix integral generating all vertices from ([TsD with equal weights 



I = J fld'^'X, d-'y exp [ - f EtrX2 _ |^^y2 + ^Nj2tT{YXaj^X, 

a=l a=l a=l 

= J d^'Fexp [ - f try2 _ |tr2 log (l ® 1 - ® ^ - ^ ® jjY)' 

(22) 

which reminds the Rostov's representation for the 0{n) matrix model in the 
dense phase |§] . We have attached the lower index to the gaussian X variable 
to weight the closed loops. Thus, 



/(^,/i) = E 



No-l,,N3 



(23) 



where J2{v} is the sum over all the diagrams. This matrix model belongs 
to the same universality class as the 0(n)-model in the dense phase. I{i^,f^) 
has a finite radius of convergence as a function of fi. Thus, the theorem is 
established. 



3 Discussion 

It is very difficult to take account of the entropy of configurations in order to 
improve the representation (p^). Therefore, let us consider the simplest self- 
consistent model, which is equivalent to the dense phase of the 0{n) matrix 
model HI: 



/o(„) = / rf^V exp [ - ytry2 - ^tr^ log (l ® 1 - ^y ® 1 - 1 ® ^v)] (24) 
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Let us suppose that this truncated model bears nevertheless some quali- 
tative features of 3-dimensional simplicial gravity. 

The critical behavior of the 0{n) matrix model is well known For 
u < 2, it describes 2-d gravity interacting with c < 1 conformal matter 

c=l-Q^^—^^; zy = -2cos7r£/o (25) 
9o 

while for ly > 2 the corresponding matter is non-critical and the model trivi- 
alizes. In this phase, the number of closed loops is proportional to the volume 
and the mean length of each remains finite. 

In the vicinity of a critical point fi^, the partition function behaves as 

/0(n) ~ (/ic - f^)^-^''' (26) 

Two main quantities of interest are 

(iV3)o(„, = f^§^ log/o(„, ^ -(2 - 7stO^ (27) 

and 



(^o)o(„) = log /o(,o ^ (2 - 7str)3^ - log 1^. - /i 



-^£(^3)o(„,+^%^log(iV3)o,„) 



(28) 



where the coefficients can be calculated explicitly using the Gaudin and Ros- 
tov's exact solution |^ 

d . 11 d'jstv 1 

— log/ic = r; -7^ = — — — n V (29) 

ou 2u + 2 ou 27igQsmTr{l — go) 

If 2 - 1/ < 1, then 

(^o)o(„) ^ \{N3)o,., + ^^7f^l°g<^3)o(„) (30) 
If > 2, then 7str = — | is independent of u and (A^o)o(„) is proportional 

to (iV3)o(„,- 

This behavior is strikingly similar to the results of numerical simulations! 
In Ref. M, a phase transition in the model (0) with respect to the fugacity 
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for the number of vertices, a, was found. In the "cold" phase, a > ac, 
the mean number of vertices, (A^o), is strictly proportional to a volume, 
{N^). In the "hot" phase, a < ac, (Nq) is a non-trivial function of (N^)- 
The analogy with the 0{n) matrix model suggests that the most probable 
scaling is| (A^o) ~ Ci{N3) + C2 log(A^3) with C2 singular at the critical point a^. 
Presumably, it is a type of this singularity that can be, in principle, calculated 
in continuum theory in order to compare predictions of both approaches. 
The obvious problem for such a hypothetical comparison is that a is a bare 
coupling. Therefore, in the lattice model, only critical points with respect to 
it could show some universal features. 

Of course, the 0{n) matrix model cannot give us precise information 
about simplicial gravity. However, it is very plausible that it has qualitatively 
the same phase structure as models (0) and (H) and may be quite instructive 
from this point of view. 

To conclude, let us express the hope that, when 3-dimensional simplicial 
gravity is solved, it will reveal as much beautiful mathematical and physical 
structure as the matrix models have been doing. 
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